Dynamical effects of the nanometer-sized polarized domains in Pb(Zni/3Nb2/3)03 
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Recent neutron scattering measurements performed on the 
relaxor ferroelectric Pb[(Zni/3Nb2/3)o.92Tio.o8]03 (PZN-8% 
PT) in its cubic phase at 500 K, have revealed an anoma- 
lous ridge of inelastic scattering centered ~ 0.2 A^^ from the 
zone center (Gehring et al, Phys. Rev. Lett. 84, 5216 (2000)). 
This ridge of scattering resembles a waterfall when plotted as 
a phonon dispersion diagram, and extends vertically from the 
transverse acoustic (TA) branch near 4 meV to the transverse 
optic (TO) branch near 9 meV. No zone center optic mode 
was found. We report new results from an extensive neutron 
scattering study of pure PZN that exhibits the same waterfall 
feature. We are able to model the dynamics of the waterfall 
using a simple coupled-mode model that assumes a strongly q- 
dependent optic mode linewidth Pi (q) that increases sharply 
near 0.2 A"^ as one approaches the zone center. This model 
was motivated by the results of Burns and Dacol in 1983, 
who observed the formation of a randomly-oriented local po- 
larization in PZN at temperatures far above its ferroelectric 
phase transition temperature. The dramatic increase in Pi 
is believed to occur when the wavelength of the optic mode 
becomes comparable to the size of the small polarized micro- 
regions (PMR) associated with this randomly-oriented local 
polarization, with the consequence that longer wavelength op- 
tic modes cannot propagate and become overdamped. Below 
Tc — 410 K, the intensity of the waterfall diminishes. At low- 
est temperatures (~ 30 K) the waterfall is absent, and we 
observe the recovery of a zone center transverse optic mode 
near 10.5 meV. 



I. INTRODUCTION 

Scientific research in the field of relaxor ferroelectrics 
has surged markedly over the past several years. Per- 
haps the primary driving force behind the rapidly in- 
creasing number of publications in this area has been 
the observation of an exceptionally large piezoelectric 
response in single crystals of the two complex per- 
ovskite systems Pb[(Mgi/3Nb2/3)i-2;Ti^]03 (PMN-xPT) 
and Pb[(Zni/3Nb2/3)i-j;Ti^]03 (PZN-xPT) for compo- 
sitions X that lie near the morphotropic phase boundary 
(MPB), which separates the rhombohedral and tetrago- 
nal regions of the phase diagram. These systems differ 
from the classic and well-studied simple perovskite ABO3 
compounds by virtue of the mixed-valence character of 



the B-site cation. In the case of PMN and PZN (a;=0), 
the i3-site is occupied by either Mg^"*" or Zn^+ , and Nb^+ 
ions with a stoichiometry of 1/3 and 2/3, respectively, 
which is required to preserve charge neutrality. This 
built-in disorder sharply breaks the translational sym- 
metry of the lattice and produces a so-called "diffuse" 
phase transition in which the dielectric permittivity e ex- 
hibits a very large and broad peak at a characteristic tem- 
perature Tmax that is also strongly frequency-dependent. 
In contrast to the normal ferroelectric parent compound 
PbTiOs (PT), in which the condensation of a transverse 
optic (TO) zone center phonon leads to a transition from 
a cubic paraelectric phase to a tetragonal long-range or- 
dered state with a non-vanishing spontaneous polariza- 
tion below a critical temperature Tc = 763 K, the proto- 
typical relaxor compound PMN exhibits a diffuse transi- 
tion at Tmax = 230 K and no spontaneous polarization 
at any temperature studied. 

Stimulated by the unusual properties of these complex 
systems, we performed a recent neutron inelastic scat- 
tering study to examine the lattice dwiamics of a high- 
quality single crystal of PZN-8%PT,El the composition 
for which th&.piezoelectric response is largest in the PZN- 
xPT system. □ The PZN-xPT phase diagram as a function 
of PT concentration x is shown in Fig. 1. The maximum 
piezoelectric activity occurs on the rhombohedral side of 
the MPB, which is indicated by the steep dashed line. 
Our measurements of the phonon dispersion of the po- 
lar TO mode in the cubic phase at 500 K revealed an 
unexpected ridge of scattering centered at a momentum 
transfer q = 0.2 A~^, measured from the zone center. 
The ridge extends in energy from ~ 4 meV to 9 meV 
such that, when plotted as a standard phonon dispersion 
diagram, it resembles a waterfall in which the TO branch 
appears to drop precipitously into the transverse acoustic 
(TA) branch. This anomalous feature, hereafter referred 
to as the waterfall, was observed in both of the two Bril- 
louin zones we surveyed (2,2,Z), and {h,h,A), which cover 
the two different symmetry directions [001] and [110]. 
Similar data taken on a large single crystal of pure PMN 
at room temperature also showed the presence of a wa- 
terfall, suggesting that this feature may be common to 
all relaxor systems.H 
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Fig. 1. Phase diagram of the PZN-PT system from Kuwata 
et al. (ref. [3]). The MPB is represented by the dashed line 
between the rhombohedral and tetragonal phases. An in- 
tervening monoclinic phase has been discovered in PZT by 
Noheda et al. (ref. [4]). 

The waterfall observed in PZN-8%PT was specu- 
lated to represent direct microscopic evidence for small, 
randomly-oriented regions of local polarization within 
the crystal that begin to condense at temperatures 
far above Tc = 450 K. The existence of these re- 
gions of local polarization was first proposed by Burns 
and Dacol in 1983|-|to describe the disorder inherent 
to relaxor systems.Q Using measurements of the op- 
tic index of refraction to study ceramic samples of 
(Pbi_3j./2Laa;)(Zrj,Tii_j,)03 £PLZT), as well as single 
crystals of PMN and PZN,eI they discovered that a 
randomly-oriented, and non-reversible local polarization 
Pd develops at a well-defined temperature Td, frequently 
referred to as the Burns temperature, several hundred 
degrees above the apparent transition temperature Tc- 
Subsequent studies have provided additional experimen- 
tal evidence supporting the existence of Td.LrEl The spa- 
tial extent of these locally polarized regions was conjec- 
tured to be of the order of several unit cells, and has 
given rise to the term "polar micro-regions," or PMR. 
Because the size of the PMR is finite, the propagation 
of long- wavelength phonons should be inhibited. This 
would then lead naturally to a heavily overdamped TO 
phonon cross section in the neighborhood of the zone 
center. Indeed, no well-defined TO phonon peaks were 
found in either PZN-8%PT or PMN for q < 0.2 A'^. 
Moreover, if one assumes the waterfall peak position in 
q gives a measure of the size of the PMR according to 
27r/(7, one obtains a value of 30 A, or about 7-8 unit 
cells, consistent with the conjecture of Burns and Dacol. 



Further evidence in support of our picture correlating 
the waterfall with the appearance of the PMR comes from 
the neutron inelastic scattering data of Naberezhnov et 
al. As shown in Fig. 2, they observed a normal TO and 
TA phonon dispersion in the closely related relaxor ferro- 
electric compound PMNj-at 800 K > = 617 K, where 
no polar regions remain.E3 Their data were obtained on 
the same crystal we studied at room temperature where 
the waterfall was observed. For PZN-8%PT, the forma- 
tion of the PMR is estimated to occur at Td ~ 790 K, 
which is well above the cubic-to-tetragonal phase transi- 
tion at Tc ~ 450 K, and also beyond the temperature at 
which the crystal begins to decompose. 
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Fig. 2. Dispersion curves of the transverse acoustic (TA) 
and the lowest-energy transverse optic (TO) modes in PMN 
measured by Naberezhnov et al. (ref. [10]) at 800 K, far above 
Ta^ 617 K. 

In this paper we present a more complete lattice dy- 
namical study of pure PZN, in which we also observe 
an anomalous ridge of scattering that closely resembles 
the waterfall seen in PZN-8%PT and PMN. PZN differs 
from PMN in that it exhibits an explicit structural phase 
transformation at 410 K from cubic to tetragonal sym- 
metry. The value of T^ax for PZN is just slightly above 
Tc = 410 |X, and thus much higher than T^ax = 230 K 
for PMN.Ilil We further present a simple mode-coupling 
model that clearly relates this unusual feature of the TO 
phonon branch (the same branch that goes soft at the 
zone center at Tc in PbTiOs) to the PMR. Our current 
study does not, however, cover the critical scattering at 
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small q below 0.075 A~^. Extensive studieSjin this area 
have alreadii-been carried out using neutront3 and x-ray 
scatteringt^EJ techniques. As we will show, this is the 
portion of the Brillouin zone where the TO mode is over- 
damped because of the PMR. 



II. EXPERIMENTAL 

Two single crystals of PZN were used in this study, 
the first, labelled PZN#1, weighs 4.2 grams and the sec- 
ond, labelled PZN#2, weighs 2.3 grams. Both crystals 
were grown using the high-temperature flux technique de- 
scribed elsewhere.til Each crystal was mounted onto an 
aluminum holder and oriented with the cubic [00 f] axis 
vertical, thereby giving access to the (HKO) scattering 
zone. They were then placed inside a vacuum furnace 
capable of reaching temperatures of up to 900 K. For 
temperatures below room temperature, a specially de- 
signed closed-cycle helium refrigerator was used to cover 
the range from 25 K to 670 K. 

All of the neutron scattering data presented here were 
obtained on the BT2 triple-axis spectrometer located at 
the NIST Center for Neutron Research. The (002) reflec- 
tion of highly-oriented pyrolytic graphite (HOPG) was 
used to monochromate and analyze the incident and scat- 
tered neutron beams. An HOPG transmission filter was 
used to eliminate higher-order neutron wavelengths. The 
majority of our data were taken holding the incident neu- 
tron energy Ei fixed at 14.7 meV (A^ = 2.36 A) while 
varying the final neutron energy Ef, and using horizontal 
beam colhmations of 60'-40'-S-40'-40' and 60'-40'-S-40'- 
open. 

Two basic types of inelastic scans were used to collect 
data on each sample. Constant-£' scans were performed 
by holding the energy transfer hu — AE — Ei~ Ef fixed 
while varying the momentum transfer Q. Constant-Q 
scans were performed by holding the momentum trans- 
fer Q = fc^ — fc/ (fc = 27r/A) fixed while varying the en- 
ergy transfer AE. Using a combination of these scans, 
the dispersions of both the TA and the lowest-energy 
TO phonon modes were mapped out at a temperature of 
500 K (still in the cubic phase, but well below the Burns 
temperature for PZN of Td ~ 750 K). 

It is worthwhile to note that the high temperature data 
were taken using neutron energy gain, i. e. by scanning 
Ef such that the energy transfer < 0. This process 
is known as phonon annihilation. At high temperatures, 
and not too large energies, the phonon thermal popula- 
tion factor makes this mode of operation feasible. More 
importantly, the instrumental energy resolution function 
contains the factor cot 26*0, where 29a is the analyzer 
(HOPG) scattering angle. Hence the energy resolution 
increases with increasing energy because 29a decreases. 
The net result is an increased neutron count rate at the 
expense of a broader linewidth, thereby making the rel- 
atively weaker transverse optic mode easier to see. 



III. INELASTIC SCATTERING DATA FROM 
PZN NEAR Q= (2,0,0) 

The waterfall, previously reported in PZN-8%Pt|3 is 
shown for PZN#1 at 500 K in Fig. 3. The data points 
represent the peak scattered neutron intensity plotted as 
a function of ?iw and g, where q = Q — G is the momentum 
transfer relative to the G = (2,0,0) Bragg reflection, 
measured along the cubic [010] symmetry direction. The 
solid dots represent data from constant-Q scans whereas 
the two open circles are data derived from constant-E' 
scans. These latter two data points signal the onset of the 
waterfall regime in which the TO phonon branch appears 
to dive into the TA branch, and are thus connected by 
a dashed line. The lowest-energy data points trace out 
the TA phonon branch along [010]. Solid lines are drawn 
through these points as a guide to the eye, and are nearly 
identical to those shown for PMN in Fig. 1. 
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Fig. 3. Solid (open) dots represent positions of peak scat- 
tered neutron intensity taken from constant-Q (constant-_E 
scans) at 500 K along the cubic [010] symmetry direction in 
PZN. Solid lines are guides to the eye indicating the TA and 
TO phonon dispersions. 



As was the case for PZN-8%PT measured at the same 
temperature, the anomalous ridge of scattering in PZN 
is centered approximately at 0.14 reciprocal lattice units 
(rlu) or 0.22 k"^ (1 rlu= 27r/a = 1.545 A"!). This is 
more clearly shown by the series of constant--E scans plot- 
ted in Fig. 4 for Tioj =0, -6, -8, and -12 meV (negative 
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values correspond to phonon annihilation as described 
above). We note that both of the -6 and -8 mcV scans 
are peaked at the same value of k {= q), thereby indicat- 
ing the presence of the waterfall. Remarkably, the meV 
cross section exhibits an enormous change in the vicin- 
ity of the waterfall g- vector. This demonstrates that the 
effects of the PMR extend in energy all the way from 
the TO branch through the TA branch and into the elas- 
tic channel. The -12 meV scan, by contrast, peaks at a 
larger momentum transfer q (~ 0.24 rlu) that lies outside 
the waterfall regime, and instead represents a genuinely 
propagating TO phonon mode. 
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Fig. 4. Constant--B scans at 0, -6, -8, and -12 meV (phonon 
annihilation) measured at 500 K on PZN crystal #1. 

It is important to realize that peaks in constant-£' 
scans do not guarantee propagating modes. This can only 
be confirmed using constant-Q scans. Therefore we show 
four representative constant-Q scans in Fig. 5 for the q 
values 0.00, 0.16, 0.28, and 0.48 rlu. For q > 0.16 rlu, two 
distinct peaks are present, corresponding to the TA and 
TO phonon modes. At 0.16 rlu, one sees the TO mode 
beginning to become damped while the intensity of the 
TA mode has become relatively much more intense. For 
q < 0.16 rlu, the TA mode was indistinguishable from the 
elastic peak. Most striking is the absence of any peak in 
the constant-Q scan for q = 0.00 rlu. There is no clear 
peak for the TO mode for q < 0.16 rlu, indicating that 



these polar optic modes are heavily overdamped. 
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Fig. 5. Constant-Q scans at 0.00, 0.24, 0.28, and 0.48 rlu 
measured at 500 K on PZN crystal #1. 



To illustrate this behavior as clearly as possible, a se- 
ries of constant-Q scans were collected in steps of 0.04 rlu 
along (2, k, 0) that spanned the entire Brillouin zone. Af- 
ter subtraction of a constant background, the data have 
been plotted as a contour plot in Fig. 6. For ease of 
comparison, the contour plot has been drawn with the 
same vertical and horizontal scales used in Fig.'s 2 and 
3. The intensity data are represented on a logarithmic 
color scale with yellow representing the highest intensity, 
and black the lowest. At large q near the zone boundary, 
the TO and TA phonon branches are readily apparent 
as broad green features. Near the waterfall regime, the 
TO branch appears to drop almost vertically into the TA 
branch with a corresponding increase in intensity. Below 
g' ~ 0.14 rlu, the intensity contours show no peak (scan- 
ning vertically) that could correspond to a propagating 
mode. The acoustic mode looks to be flat below 4 meV 
only because the color scale covers a limited intensity 
range above which all data arc colored yellow. The inten- 
sity range was limited on purpose to make the waterfall 
feature more visible. Even so, one can notice a very rapid 
increase in the TA phonon intensity in the vicinity of the 
waterfall, suggesting that a transfer of intensity might be 
occurring between the TO and TA modes. To examine 
this possibility, we performed model calculations using a 
simple mode-coupling scheme to attempt to simulate the 
anomalous features reported here in the PZN compound. 
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Fig. 6. Contour map of the background subtracted scatter- 
ing intensity from PZN at 500 K measured near (200). The 
intensity is indicated by a logarithmic color scale that is lim- 
ited to a narrow range in order to better show the waterfall. 
Yellow is the most intense. 



IV. MODE-COUPLING MODEL 

During the course of our experiments on relaxor fer- 
roelectrics, a simple but effective model describing the 
anomalous low-frequency lattice dynamics was developed 
assuming a coupling between the TO and TA phonon 
modes. For the case of neutron energy loss, the scat- 
tering intensity distribution / for two interacting modes 
with frequencies fli arui.rJ2, and widths Fi and F2, is 
given by the expressionEj 



[n{uj) + l]- 



ml-co^)B-T,A]Fi), 



IXB F1F2 



where A and B are given by 



A 
B 



Ti{nl-oj'^) + T2{nl 



(1) 



(2) 



and n(w) is simply the Bose factor \e^^/^BT) _ 
The quantities Fi^2 are the structure factors of modes 
1 and 2, and A is the coupling strength between the two 



modes. Extensive model calculations based on this equa- 
tion, which has been shown to describe the-lpehavior of 
coupled-phonon cross sections quite welljl30 have been 
performed in collaboration with K. Ohwada at the Insti- 
tute of Solid State Physics, University of Tokyo.ta 
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Fig. 7. Model simulations assuming a coupled-mode in- 
tensity distribution and a strongly g-dependent TO phonon 
linewidth Fi. Three constant-Q scans are shown correspond- 
ing to q = 0.20, 0.15, and 0.10 rlu, with q = 0.15 rlu taken to 
be the reciprocal space position of the anomalous waterfall 
feature. 



The essential physics behind the mode-coupled model 
description of the waterfall is built into the linewidth of 
the TO mode Fi, which is assumed to become sharply 
q-dependent as the polar micro-regions begin to form 
at the Burn's temperature Td- If we suppose that the 
PMR have an average diameter given by 27r/q^/, where 
qwf represents the reciprocal space position of the wa- 
terfall, then those optic phonons having q < q^f will 
not be able to propagate easily because their wavelength 
exceeds the average size of the PMR. These polar lat- 
tice vibrations are effectively impeded by the boundary 
of the PMR. This idea has also been discussed by Tsu- 
rumi et a/O The simplest way to simulate this situation 
is to assume a sudden and steep increase in Fi(g) at 
For this purpose, a Fermi-distribution function olq works 
extremely well. Fig. 7 shows several model constant-Q 
simulations based on this assumption, using the values 
g^/ = 0.15 rlu, A = 10, Fi = 1, F2 = 4, and F2 = 1. For 



5 



simplicity and purposes of illustration, the dispersions of 
both optic and acoustic modes were ignored by holding 
the parameters fii,2 fixed at 10 and 5 meV, respectively, 
over the entire Brillouin zone. Similarly, instrumental 
resolution effects were not included. 

For q > qyjf one observes two broad peaks, as expected. 
At momentum transfers q < q-wf, however, the optic 
mode becomes highly overdamped and its profile extends 
in energy below that of the acoustic mode. Alongside 
each constant-Q scan is shown the corresponding value 
of Fi used in the simulation. The waterfall thus rep- 
resents the crossover between a high-g regime, in which 
one observes two well-defined peaks corresponding to two 
propagating modes, and a low-g regime, in which one ob- 
serves an overdamped optic mode plus an acoustic peak. 
This simple model cross section describes all of the exper- 
imental observations very well. Indeed, one can favorably 
compare the simulated scan at g = 0.15 rlu in Fig. 7 with 
the experimental scan for Q = (2, 0.16, 0) shown in Fig. 5. 
One can see now that the waterfall is not a TO phonon 
dispersion at all. Instead, it is simply a redistribution of 
the optic mode profile that is caused by the PMR which 
force a sudden change in the optic mode linewidth at a 
specific q related to the average size of the PMR. 
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Fig. 8. Model simulations of two constant-iJ scans at and 
7.5 meV. The g-dependence of the TO linewidth Fi is shown 
in the top panel for ease of comparison. 



scattering cross section within the coupled-mode model, 
the g-dependence of the optic mode linewidth Fi is shown 
in Fig. 8 along with two simulated constant- scans at 
and 7.5 meV. It is apparent that the sharp increase in Fi 
has a corresponding and pronounced effect on both cross 
sections in the vicinity of q^f- The jump in the elastic 
cross section occurs because of the sum rule for the scat- 
tering. Hence the extreme broadening of the TO mode 
that occurs near q^f extends underneath the TA branch 
and into the elastic channel. We emphasize that despite 
the fact that no resolution effects, Bragg or background 
terms were included in these calculations, the agreement 
with experiment is surprisingly good. Indeed, the cross 
section calculated for meV is striking in its similarity to 
the experimentally measured elastic cross section shown 
in the top panel of Fig. 4, and thus lends strong credence 
to the coupled-mode model. 



V. TEMPERATURE DEPENDENCE 

Our data on the temperature dependence of the wa- 
terfall in PZN are limited by the fact that PZN decom- 
poses at a temperature less than the Burns temperature. 
This prevents us from performing measurements above 
Td which are of interest because we wish to monitor the 
recovery of zone center and low-g TO phonons at temper- 
atures where the polarized micro-regions no longer exist. 
PMN, by contrast, has a lower Td- Hence Naberezh- 
nov et al. were able to measure the TO phonon disper- 
sion at 800 K >-,Td where they observed a normal TO 
phonon branchJlj Selected phonons were briefly studied 
at 860 K in PZN-8%PT.EI There the scattering intensity 
associated with the waterfall increased (as expected from 
the Bose temperature factor) beyond that measured at 
500 K, and the peak position of the waterfall, measured 
using constant--E scans, shifted in towards the zone cen- 
ter, i. e. from (2,0.14,0) to (2,0.10,0). Constant-Q data at 
(2,0.10,0) show a broad peak around 10 meV. These data 
suggest that around 860 K the PMR have already begun 
to form inside of a "uniform," or unpolarized crystal. A 
complete understanding of the behavior in this tempera- 
ture range must await the further study of PMN. 

Below the PZN tetragonal-to-rhombohedral phase 
transition at 410 K the scattering intensity associated 
with the waterfall decreases gradually, and is almost com- 
pletely gone at 100 K. This behavior is summarized in 
the inset to Fig. 9 where the raw intensity (no back- 
ground subtraction) at Q = (2, —.12, 0) and Ticu = 7 meV 
is plotted versus temperature. Figure 9 itself shows a 
constant-Q scan at (2,0,0) and 25 K which clearly shows 
the presence of a zone center TO mode. 



In order to illustrate the basic characteristics of the 
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Fig. 9. Constant-Q scan at (2,0,0) and 25 K taken on PZN 
crystal #2 showing the recovery of the TO zone center mode 
at low temperature. Inset gives the temperature dependence 
of the raw scattering intensity measured at Q = (2,-0.12,0) 
at 7 meV (no background subtraction). 



at low q. A brief review of the waterfall and related neu- 
tron inelastic scattering measurements on various relaxor 
systems has recently been submitted for publicationE^I It 
is hoped that our results with stimulate theoretical cal- 
culations of the precise q-depcndcnce of Fi to allow for 
comparison with experiment. It would be very interest- 
ing to investigate the effects of an applied electric field 
E on the waterfall, and we have already begun measure- 
ments to this end. 
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VI. DISCUSSION 

We have discovered an anomalous ridge of scattering 
positioned roug hly 0.2 from the zone center in the 
relaxor systems PZN-8%PT, PZN, and also PMN in their 
respective cubic phases. This scattering was observed in 
various different Brillouin zones including (220) and (004) 
in PZN-8%PT, and (200), (300), and (220) in PZN. The 
ridge of scattering is extended in energy in such a way as 
to make the lowest-lying TO phonon branch appear to 
drop precipitously into the TA phonon branch at a spe- 
cific value of q, thereby resembling a waterfall. Equally 
remarkable is the manifestation of a huge elastic cross 
section near the same value of q. An unusual feature 
has also been observed by Lushnikov et al. for pure PMN 
using Raman and neutron inelastic scattering techniques 
at 77 K over an extended energy range from 5.2 meV 
to 6.8 meV, which nearly coipcides with that of the wa- 
terfall feature reported here.Ej This was interpreted as 
evidence of a fracton, i. e. an excitation on a fractal lat- 
tice. We have shown that a simple coupled-mode model 
can be used effectively to describe our neutron inelastic 
data by postulating a sharply g-dependent TO phonon 
linewidth Fi that exhibits a large step-like increase with 
decreasing q upon reaching the waterfall wavevector qyjf. 
In this way we are able to relate the waterfall directly to 
the presence of nanometer-sized polarized micro-regions 
in the crystal which serve to damp the polar TO modes 
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